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Abstract 

Group A rotavirus (RVA) infections form a major public health problem, especially in low-income countries like the 
Democratic Republic of the Congo (COD). However, limited data on RVA diversity is available from sub-Saharan Africa in 
general and the COD in particular. Therefore, the first aim of this study was to determine the genetic diversity of 99 RVAs 
detected during 2007-2010 in Kisangani, COD. The predominant G-type was G1 (39%) and the most predominant P-type 
was P[6] (53%). A total of eight different G/P-combinations were found: G1P[8] (28%), G8P[6] (26%), G2P[4] (14%), G12P[6] 
(13%), G1P[6] (11%), G9P[8] (4%), G4P[6] (2%) and G8P[4] (1%). The second aim of this study was to gain insight into the 
diversity of P[6] RVA strains in the COD. Therefore, we selected five P[6] RVA strains in combination with the G1, G4, G8 (2x) 
or G12 genotype for complete genome analysis. Complete genome analysis showed that the genetic background of the 
G1P[6] and G12P[6] strains was entirely composed of genotype 1 (Wa-like), while the segments of the two G8P[6] strains 
were identified as genotype 2 (DS-1-like). Interestingly, all four strains possessed a NSP4 gene of animal origin. The analyzed 
G4P[6] RVA strain was found to possess the unusual G4-P[6]-I1-R1-C1-M1-A1-N1-T7-E1-H1 constellation. Although the 
majority of its genes (if not all), were presumably of porcine origin, this strain was able to cause gastro-enteritis in humans. 
The high prevalence of unusual RVA strains in the COD highlights the need for continued surveillance of RVA diversity in the 
COD. These results also underline the importance of complete genetic characterization of RVA strains and indicate that 
reassortments and interspecies transmission among human and animal RVAs strains occur regularly. Based on these data, 
RVA vaccines will be challenged with a wide variety of different RVA strain types in the COD. 
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Introduction 

Group A rotaviruses (RVAs) are the world's leading cause of 
severe diarrhea in children <5 years of age and form a major 
public health problem, with an estimated 453,000 deaths per year 
worldwide. More than 50% of these deaths occur in Africa 
(232,000) and 7% of all RVA related deaths occur in the 
Democratic Republic of the Congo (COD) (32,653; 95% CI 
27,804-37,699) [1]. In fact, the COD, together with India, 
Nigeria, Pakistan and Ethiopia, is one of the countries that bear 
the highest RVA mortality rates worldwide. 

Despite this high burden of disease, only limited regional and 
country specific data on RVA genotypes diversity is available from 
sub-Saharan Africa in general and the COD in particular. The last 
report of the African Rotavirus Surveillance Network reported the 
distribution of predominant RVA genotypes among hospitalized 



children aged <5 years between 2007 and 201 1 from 16 African 
countries [2]. The major G- and P-genotype combinations 
detected in this study were G1P[8] (18.4%), G9P[8] (11.7%), 
G2P[4] (8.6%), G2P[6] (6.2%), G1P[6] (4.9%), G3P[6] (4.3%), 
G8P[6] (3.8%) and G12P[8] (3.1%). Worldwide the majority of 
RVA strains infecting humans belong to the G 1 -G4, G9 and G 1 2 
G-genotypes, in combination with the P[4] (for G2) and P[8] (for 
Gl, G3, G4, G9 and G12) genotypes. When focusing on the 
differences in the genetic setup of RVA strains circulating in 
Africa, compared to the ones circulating in industrialized 
countries, multiple studies have reported the high detection rate 
of RVA strains containing the P[6] and G8 genotypes. Both 
genotypes are believed to be of animal origin: P[6] RVA strains 
(most likely) originated from pigs, while G8 is a genotype often 
detected in cows or other members of the mammalian order of the 
Artiodactyla [3-10]. The majority of known human G8 RVAs 
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strains have been described all over the African continent in 
combination with a large number of VP4 specificities and with 
different genetic backgrounds [11-23]. 

Two human RVA vaccines, Rotarix (GlaxoSmithKline) and 
RotaTeq (Merck) have been licensed and proven safe in large 
clinical trials in the Americas and Europe [24,25]. Data from 
clinical trials conducted in five African countries have shown a 
lower vaccine immunogenicity and efficacy compared to the ones 
observed in America and Europe [26-28]. Despite this observed 
lower vaccine efficacy, vaccine introduction in Africa will prevent 
a greater number of RVA related deaths due to the substantially 
larger mortality rates in Africa compared to the Americas and 
Europe. Between August 2009 and May 2014, 21 African 
countries (in chronological order: South-Africa, Botswana, Mo- 
rocco, Sudan, Ghana, Rwanda, Malawi, Tanzania, The Gambia, 
Burkina Faso, Ethiopia, Libya, Zambia, Burundi, Mali, Camer- 
oon, Sierra Leone, The Republic of the Congo, Angola, 
Madagascar and Zimbabwe) have implemented RVA vaccination 
into their national immunization program [29,30]. 

What factors determine the lower vaccine efficacy in Africa is 
not yet completely understood but both host (malnutrition, 
competing enteric pathogens, differences in histo-blood group 
antigens, the level of anti-rotavirus antibodies in breast milk and a 
lower immunological response) and viral (a higher RVA strain 
diversity) characteristics have been proposed to (partially) explain 
this observation [31-33]. 

Specific data on RVA strain diversity from the Democratic 
Republic of the Congo are scarce with only two currently available 
studies dealing with the RVA strain diversity in the COD. The 
first study analyzed the complete genomes of two G8 strains 
(RVA/Human-wt/COD/DRC86/2003/G8P[6] and RVA/Hu- 
man-wt/COD/DRC88/2003/G8P[8]). These strains both pos- 
sess a DS-l-like genotype constellation and are the only complete 
genomes of Congolese RVA strains available today [13]. The 
second study described the epidemiology of RVA associated 
gastroenteritis in Kinshasa, COD, from July to October in 2003, 
2004, and 2005. The predominant genotypes in 2003 were 
G8P[6] and G8P[8] strains while in 2004 and 2005 G1P[6] strains 
were predominant [14]. 

The number of studies conducted in sub-Saharan Africa 
involving the analyses of complete genomes of human RVA 
strains isolated in Africa is also limited. To our knowledge the 
complete genotype constellation of 87 human RVA strains isolated 
in sub-Saharan Africa are retrievable from Genbank: 28 of the 87 
RVA strains were isolated in Malawi, 22 in Cameroon, 20 in 
South-Africa, 4 in Kenya, 3 in Nigeria, 3 in Cote d'lvoire, 3 in 
Zambia, 2 in the COD, 1 in Zimbabwe and 1 in Ethiopia 
[11,13,34—38]. Compared to the hundreds of complete genomes 
characterized from RVAs isolated in high-income countries, and 
considering the large number of co-circulating G/P genotype 
combinations, this number is rather low. 

The present study aimed to determine the genetic diversity of 
RVAs detected during 2007-2010 from Kisangani, COD. In 
addition we selected five P[6] RVA strains in combination with the 
different observed G-genotypes (Gl, G4, 2xG8 and G12) for 
complete genome analyses to gain insight into the diversity, origin 
and evolution of P[6] RVA strains in the COD. 

Methods 

Stool sample collection and detection of RVA 

Fecal specimens were collected between May 2007 and 
December 2010 from children <5 years old admitted with 
gastroenteritis (GE) to the pediatric wards of three hospitals in 



Kisangani (University Hospital, General Referral Hospital Maki- 
sao of Kisangani - both public hospitals - and Village de Pediatric, 
a private hospital center) and screened by an immunochromato- 
graphic antigen test (Rota-CIT, BioConcept, Belgium) for the 
presence of RVA antigen. Fecal samples containing RVA were 
collected and shipped to the Rega Institute for Medical Research 
using chromatography paper strips as previously described [39] . 

Rotavirus genotyping 

The paper strips were inserted into an Eppendorf tube with 
500 ul of universal transport medium (Copan Diagnostics, 
Corona, CA, USA) and thoroughly squeezed using sterile forceps. 
An aliquot of 140 |il of the squeezed eluate was used for RNA 
extraction using the QIAamp Viral RNA mini kit (QIAGEN/ 
Westburg, Leusden, The Netherlands) according to the manufac- 
turer's instructions. The G- and P-genotypes of 99 antigen-positive 
samples were characterized by reverse-transcription polymerase 
chain reaction (RT-PCR) using the QIAGEN OneStep RT-PCR 
kit (QIAGEN/Westburg) using primers Beg9 and End9 (VP7), 
and VP4_1-17F and Con2 (VP4). Primer sequences are shown in 
table S 1 . PCR products were run on a polyacrylamide gel, stained 
with EtBr and visualized under UV-light. The positive samples 
were purified with ExoSAP-IT (Affymetrix, Santa Clara, CA, 
USA), and Sanger sequenced using the forward primers (Beg9 or 
VP4_1-17F) with the ABI PRISM BigDye Terminator cycle 
sequencing reaction kit (Applied Biosystems, Foster City, CA, 
USA). 

Sequencing complete genomes 

For the complete genomes, RT-PCR was carried out with an 
initial reverse transcription step at 50°C for 30 min followed by a 
PCR activation step at 95°C for 15 min, 35 cycles of amplification, 
and a final extension step for 10 min at 70°C in a Biometra T3000 
thermocycler (Biometra, Westburg BV, Netherlands). For the gene 
segments encoding VP6, VP7, NSP1, NSP2, NSP3, NSP4 and 
NSP5, the amplification cycle conditions were as follows: 30 s at 
94°C, 30 s at 45°C, and 2 min at 72°C. For the larger segments 
encoding VP1, VP2, VP3 and VP4 the cycle conditions were 30 s 
at 94°C, 30 s at 48°C, and 6 min at 72°C. Primers used to amplify 
the VP1, VP2, VP3, VP6, NSP1, NSP2, NSP3, NSP4 and NSP5 
gene segments are described in table S 1 . Primer walking was used 
for full-length sequencing of the longest genes (VP1-VP4). The 5' 
and 3 ' terminal sequences of the 11 genomic segments were 
determined for all five strains using a modified RACE technique as 
previously described [40]. 

Sequence alignment and phylogenetic analysis 

The Mega 5.10 software was used for phylogenetic and 
molecular evolutionary analyses [41]. To compute evolutionary 
distances we used the p-distance nucleotide substitution model. In 
phylogenetic analyses we opted for maximum likelihood phyloge- 
netic analyses using the complete deletion function in order to 
include partially sequenced genes. Bootstrap resampling analysis 
(500 replicates) was performed to measure the reliability of the tree 
topologies. 

Nucleotide sequence accession numbers 

Accession numbers from sequences of this study: KJ870690- 
KJ870932. 

Ethics Statement 

Pending the installation of an ethical committee in the Oriental 
Province of the COD and at the Kisangani hospitals, the study was 
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Figure 1. Numbers of children hospitalized with GE (blue bars) and children hospitalized with GE tested RVA positive (red bars) per 
month from May 2007 until December 2010. 

doi:10.1371/journal.pone.0100953.g001 



assessed and approved by the Provincial Health Officer (the 
highest ranked medical authority in the Oriental Province) and the 
Director's Board of UH Kisangani, General Reference Hospital of 
Makiso and Village de Pediatrie, respectively. All samples were 
collected in a standard of care setting for therapeutic or diagnostic 
reasons. There has been no influence on the treatment of the 
patients based on the results of the RVA detection test, and 
analyses for this study were performed after installment of the 
patients' treatment. Data were fully anonymized at collection. 



Results 

1. Rotavirus incidence from May 2007 to December 2010 
RVA detection rate and study data. A total of 524 fecal 
specimens were collected over the entire study period of which 154 
samples (29%) tested positive with the immunochromatographic 
antigen test. Figure 1 shows the number of children hospitalized 
with gastroenteritis (GE) and the number of RVA positive cases 
over the entire study period. There were clear fluctuations in the 
number of RVA positive cases and in the number of hospitalized 




Figure 2. Number of hospitalized GE cases (blue bars) and number of hospitalized GE cases tested RVA positive (red bars) pooled 
per month over the entire study period. The green line indicates the percentage RVA positive cases tested in each month. 
doi:1 0.1 371 /journal. pone.0100953.g002 
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Figure 3. Geographic location of the Democratic Republic of 
the Congo on the African continent; location of the city of 
Kisangani in the COD together with the genotype distribution 
of the 99 analyzed RVA strains. Percentages of the different 
genotypes are indicated on the left. The number of strains per 
genotype is shown on the right hand side. 
doi:10.1371/journal.pone.0100953.g003 

GE cases per month with a three-month long period (April^June 
2008) during which no samples were collected because of stock 
rupture. Duration of hospitalizations associated with RVA varied 
between 1 and 1 0 days and occurred among children aged 1 week- 
2 years, with a mean age of 10 months. 

To determine a possible seasonal variability in RVA disease in 
Kisangani, all data shown in figure 1 were pooled per month and 
plotted together with the percentage of RVA positive cases in 
figure 2. The proportion of RVA positive cases among children 
hospitalized with GE ranged between 18.8% and 37.1%, with an 
average of 29.0%. RVA infections occurred year-round with three 
periods of increased RVA detection rates: February, April-May 
and September-November. 

2. Genotype distribution 

G- and P- genotyping of rotavirus strains. Ninety-nine 
samples out of 154 samples were selected, based on the availability 
of sufficient stool, for further G and P-genotyping. In the analyzed 
samples the predominant G-type was Gl (detected in 39% of 
specimens) and the most predominant P-type was P[6] (53%). A 
total of eight different G/P-combinations were found: G1P[8] 
(N = 28), G8P[6] (N = 26), G2P[4] (N=14), G12P[6] (N=13), 
G1P[6] (N=ll), G9P[8] (N = 4), G4P[6] (N = 2) and G8P[4] 
(N= 1) (figure 3). Only 46% of the isolates were found with the 
most typical human RVA genotype combinations (G1P[8], 
G2P[4] and G9P[8]). In contrast, 53% of the strains were found 
with the P[6] genotype, in combination with the VP7 genotypes 
G8, Gl, G12 or G4. 

To investigate the genetic diversity within the VP4 and VP7 
genes isolated in this study phylogenetic trees based on the partial 
nucleotide sequences of both genes were constructed (figures 4—5). 
For VP4, the majority (53%) of the characterized RVA strains 
clustered within the P[6] genotype, more specifically in two 
previously defined lineages, P[6]-I and P[6]-V, both of them 
containing human and porcine strains [10]. The Congolese RVA 
strains belonging to P[6]-I lineage analyzed in this study can be 
divided in three different clusters, two large and one smaller cluster 
(figure 4). The first large cluster contains the majority (24/26) of 
G8P[6] strains while the second big cluster contains all G12P[6] 
strains. G1P[6] strains were found in both the G8P[6] as the 



G12P[6] clusters (figure 4). The third, smaller cluster contains only 
2 G8P[6] strains: RVA/Human-wt/COD/KisB101/2007/ 
G8P[6] and RVA/Human-wt/COD/KisB605/2008/G8P[6], 
the only two G8P[6] strains isolated before 2009. Both G4P[6] 
strains detected in this study (KisB332 and KisB503) belonged to 
the P[6]-V lineage, a lineage containing mostly porcine or human 
strains believed to be of (partial) porcine RVA origin. Both samples 
were isolated in 2008, only one month apart from each other in 
the same hospital and shared 99.7% nucleotide (nt) identity in 
their VP4 genes. 

All P[4] strains were found in combination with the G2 
genotype, except for a single G8 strain (KisB566). All G2P[4] 
strains were very closely related to each other despite the fact that 
they were isolated over the entire study period. The fact that the 
VP4 gene of KisB566 was almost identical to the G2P[4] strains 
suggests a rather recent reassortment event involving the VP7 
gene. The P[8] strains from Kisangani showed limited diversity, 
with almost identical G1P[8] strains circulating in the Kisangani 
area between 2007 and 2009. The only characterized G1P[8] 
strain isolated in 2010 (RVA/Human-wt/COD/KisB563/2010/ 
G1P[8]) was found to cluster closer to the four detected G9P[8] 
strains than to the other G1P[8] strains characterized in this study, 
again suggesting a recent reassortment event. 

For VP7, we found limited diversity within each G-genotype, 
especially within G2, G4, G9 and G12 were we found at least 
97.9%, 99.2%, 98.7% and 98.2% of nt similarity, respectively. 
Within the Gl genotype, two clusters could be defined, one of 
which contains only a single Kisangani strain: RVA/Human-wt/ 
COD/KisB563/2010/GlP[8] (figure 5). This strain is the only 
G1P[8] RVA strain isolated in 2010, and shows 94.7% -95.7% nt 
similarity with G1P[8] strains isolated between 2007-2009. The 
large Gl cluster showed strains with P[6] and P[8] intermingled, 
further indicating the frequent occurrence of reassortments. A 
clustering according to year of isolation can also be found within 
the G8 RVA strains. More specifically, strains isolated in 2007 or 
2008 and strains isolated in 2009 or 2010 differ 3.6% to 4.7% on 
the nucleotide level. 

Both for VP4 and VP7, all strains isolated in Kisangani were 
closely related to RVAs present in the global RVA data collection, 
in particular to strains isolated from other African patients. For 
more detail, extended trees for both VP4 and VP7 were added as 
supporting information (Figure SI for VP4 and Figure S2 for 
VP7). Both trees contain twelve representative strains from this 
study together with RVA VP4 or VP7 sequences available in 
Genbank. 

3. Complete genome analysis 

Genotype constellations of whole genomes. To gain 
insight into the diversity, origin and evolution of P[6] RVA 
strains in the COD we selected P[6] strains, based on the 
phylogenetic clustering in the VP7 tree and on the viral load of the 
samples. In total, the nucleotide sequences of the eleven genes of 
five P[6] RVA strains with the different G-genotypes detected hi 
this study (Gl, G4, 2xG8 and G12) were completely determined, 
except from the first 33 nucleotides of the VP2 gene of strain 
KisB332 which could not be completed due to the lack of sample 
available to complete the analysis. Complete genotype assignments 
of the five fully sequenced P[6] strains together with a number of 
other representative strains are shown in table 1 . 

The genetic background of the G1P[6] and G12P[6] strains was 
entirely composed of genotype 1 (Wa-like). To our knowledge, no 
strain with the same genotype constellation as RVA/Human-wt/ 
COD/KisB504/2009/GlP[6] has previously been reported from 
Africa. However, the same constellation was reported in a study 
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RVA/Human-wt/COD/KisB538/2009/G8P[6] 
RVA/Human-wt/COD/KisB543/2009/G8P[6] 
RVA/Human-wt/COD/KisB537/2009/G1P[6] 
RVA/Human-wt/COD/KisB535/2010/G8P[6] 
RVA/Human-wt/COD/KisB532/2010/G8P[6] 
RVA/Human-wt/COD/KisB530/2010/G8P[6] 
RVA/Human-wt/COD/KisB525/2009/G8P[6] 
RVA/Human-wt/COD/KisB524/2009/G8P[6] 
RVA/Human-wt/COD/KisB544/2009/G8P[6] 
RVA/Human-wt/COD/KisB545/2009/G8P[6] 
RVA/Human-wt/COD/KisB549/2010/G8P[6] 
RVA/Human-wt/COD/KisB551/2010/G8P[6] 
99 . RVA/Human-wt/COD/KisB567/2010/G8P[6] 
RVA/Human-wt/COD/KisB568/2010/G8P[6] 
RVA/Human-wt/COD/KisB529/2010/G8P[6] 
RVA/Human-wt/COD/KisB559/2010/G8P[6] 
RVA/Human-wt/COD/KisB560/2010/G8P[6] 
RVA/Human-wt/COD/KisB561/2010/G8P[6] 
RVA/Human-wt/COD/KisB562/2010/G8PI6] 
RVA/Human-wt/COD/KisB564/2010/G8P[6] 
ARVA/Human-wt/COD/KisB565/2010/G8P[6] 
, RVA/Human-wt/COD/KisB550/2010/G1P[6] 
lRVA/Human-wt/COD/KisB557/2010/G1P[6] 
RVA/Human-wt/COD/KisB555/2010/G8P[6] 
,RVA/Human-wt/COD/KisB540/2010/G8P[6] 
I RVA/Human-wt/COD/KisB553/2010/G8P[6] 
lARVA/Human-wt/COD/KisB554/2010/G8P[6] 

,RVA/Human-wt/COD/KisB101/2007/G8P[6] 

~57l-RVA/Human-wt/COD/KisB605/2008/G8P[6] 
RVA/Human-wt/COD/KisB516/2008/G1P[6] 
ARVA/Human-wt/COD/KisB521/2008/G12P[6] 
RVA/Human-wt/COD/KisB328/2009/G12P[6] 
■ RVA/Human-wt/COD/KisB508/2009/G12P[6] 
, RVA/Human-wt/COD/KisB336/2008/G1 P[6] 
lRVA/Human-wt/COD/KisB609/2008/G12P[6] 
lRVA/Human-wt/COD/KisB323/2008/G12P[6] 
RVA/Human-wt/COD/KisB511/2009/G12P[6] 
-RVA/Human-wt/COD/KisB539/2009/G12P[6] 
RVA/Human-wt/COD/KisB610/2008/G12P[6] 
tRVA/Human-wt/COD/KisB608/2008/G12P[6] 
lRVA/Human-wt/COD/KisB613/2008/G12P[6] 
|RVA/Human-wt/COD/KisB607/2008/G12P[6] 
rRVA/Human-wt/COD/KisB526/2009/G12P[6] 
7 JlRVA/Human-wt/COD/KisB534/2009/G12P[6] 
-RVA/Human-wt/COD/KisB507/2009/G1P[6] 
RVA/Human-wt/COD/KisB512/2009/G1P[6] 
RVA/Human-wt/COD/KisB405/2009/G1P[6] 
RVA/Human-wt/COD/KisB501/2009/G1P[6] 
▲ RVA/Human-wt/COD/KisB504/2009/G1 P[6] 
RVA/Human-wt/COD/KisB509/2009/G1P[6] 



l RVA/Human-wt/COD/KisB332/2008/G4P[6] 



i^RVA/Human-wt/COD/KisB603/2008/G4P[6] 



P[6] 



P[6]-l 



P[6]-V 



0.05 nt subst./site 



RVA/Human-wt/COD/KisB547/2009/G2P[4] 
RVA/Human-wt/COD/KisB558/2009/G2P[4] 
RVA/Human-wt/COD/KisB542/2009/G2P[4] 
RVA/Human-wt/COD/KisB541/2009/G2P[4] 
RVA/Human-wt/COD/KisB523/2009/G2P[4l 
RVA/Human-wt/COD/KisB566/2010/G8P[4] 
RVA/Human-wt/COD/KisB548/2009/G2P[4] 
RVA/Human-wt/COD/KisB536/2009/G2P[4] 
RVA/Human-wt/COD/KisB527/2009/G2P[4] 
— RVA/Human-wt/COD/KisB531/2009/G2P[4] 
|-RVA/Human-wt/COD/KisB325/2009/G2P[4] 
-RVA/Human-wt/COD/KisB604/2008/G2P[4] 
RVA/Human-wt/COD/KisB102/2007/G2P[4] 
-RVA/Human-wt/COD/KisB107/2007/G2P[4] 
-RVA/Human-wt/COD/KisB546/2009/G2P[4] 
RVA/Human-wt/COD/KisB563/201 0/G1 P[8] 
RVA/Human-wt/COD/KisB502/2009/G9P[8] 
RVA/Human-wt/COD/KisB556/2009/G9P[8] 
RVA/Human-wt/COD/KisB331/2009/G9P[8] 
RVA/Human-wt/COD/KisB533/2010/G9P[8] 
,— RVA/Human-wt/COD/KisB515/2009/G1P[8] 
r RVA/Human-wt/COD/KisB612/2008/G1 P[8] 
RVA/Human-wt/COD/KisB103/2007/G1P[8] 
■ RVA/Human-wt/COD/KisB401/2009/G1 P[8] 
RVA/Human-wt/COD/KisB1 06/2007/G1 P[8] 
RVA/Human-wt/COD/KisB104/2007/G1P[8] 
RVA/Human-wt/COD/KisB333/2009/G1P[8] 
RVA/Human-wt/COD/KisB403/2009/G1P[8] 
■ RVA/Humarvwt/COD/KisB505/2009/G1 P[8] 
RVA/Human-wt/COD/KisB327/2007/G1P[8] 
RVA/Human-wt/COD/KisB503/2009/G1P[8] 
RVA/Human-wt/COD/KisB611/2008/G1P[8] 
RVA/Human-wt/COD/KisB601/2008/G1P[8] 
RVA/Human-wt/COD/KisB616/2008/G1P[8] 
RVA/Human-wt/COD/KisB615/2008/G1P[8] 
RVA/Human-wt/COD/KisB614/2008/G1P[8] 

RVA/Human-wt/COD/KisB510/2009/G1P[8] 
RVA/Human-wt/COD/KisB324/2008/G1P[8] 
RVA/Human-wt/COD/KisB404/2009/G1P[8] 
RVA/Human-wt/COD/KisB506/2009/G1P[8] 
RVA/Human-wt/COD/KisB513/2008/G1P[8] 
RVA/Human-wt/COD/KisB514/2008/G1P[8] 
RVA/Human-wt/COD/KisB518/2008/G1P[8] 
RVA/Human-wt/COD/KisB519/2008/G1P[8] 
RVA/Human-wt/COD/KisB520/2008/G1P[8] 
RVA/Human-wt/COD/KisB522/2008/G1P[8l 
RVA/Hlimai>wt/COD/KisB602/2008/G1P[8] 



P[4] 



71 :, 



P[8] 



Figure 4. Maximum likelihood tree based on 496 nucleotides of the VP4 encoding gene of the 99 RVA strains characterized in this 
study. Tree was rooted using RVA/Human-tc/JPN/AU-1/1982/G3P[9] as an outgroup. Bootstrap values (500 replicates) above 70 are shown. The 
strains completely analyzed in this study are marked with a black triangle. 
doi:10.1371/journal.pone.0100953.g004 
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RVA/Human-wt/COD/KisB614/2008/G1P[8] 
RVA/Human-wt/COD/KisBB15/2008/G1P[8] 
RVA/Human-wt/COD/KisB602/2008/G1P[8] 
RVA/Human-wt/COD/KisB557/2010/G1P[6] 
RVA/Human-wt/COD/KisB550/2010/G1P[6] 
RVA/Human-wt/COD/KisB537/2009/G1P[6] 
RVA/Human-wt/COD/KisB522/2008/G1P[8] 
RVA/Human-wt/COD/KisB520/2008/G1P[8] 
RVA/Human-wt/COD/KisB519/2008/G1P[8] 
RVA/Human-wt/COD/KisB518/2008/G1P[8] 
RVA/Human-wt/COD/KisB516/2008/G1P[6] 
RVA/Human-wt/COD/KisB510/2009/G1P[8] 
RVA/Human-wt/COD/KisB506/2009/G1P[8] 
RVA/Human-wt/COD/KisB505/2009/G1P[8] 
RVA/Human-wt/COD/KisB503/2009/G1P[8] 
RVA/Human-wt/COD/KisB401/2009/G1P(8] 
RVA/Human-wt/COD/KisB336/2008/G1P[6] 
RVA/Human-wt/COD/KisB333/2009/G1P[8] 
RVA/Human-wt/COD/KisB327/2007/G1P[8] 
RVA/Human-wt/COD/KisB324/2008/G1P[8] 
r RVA/Human-wt/COD/KisB61 1/2008/G1 P[8] 
RVA/Human-wt/COD/KisB601/2008/G1P[8] 
lRVA/Human-wt/COD/KisB616/2008/G1P[8] 
RVA/Human-wt/COD/KisB514/2008/G1P[8] 
RVA/Human-wt/COD/KisB513/2008/G1P[8] 
RVA/Human-wt/COD/KisB403/2009/G1P[8] 
RVA/Human-wt/COD/KisB404/2009/G1P[8] 
• RVA/Human-wt/COD/KisB515/2009/G1P[8] 
RVA/Human-wt/COD/KisBB12/2008/G1P[8] 
RVA/Human-wt/COD/KisB1 04/2007/G1 P[8] 
RVA/Human-wt/COD/KisB106/2007/G1P[8] 
RVA/Human-wt/COD/KisB1 03/2007/G1 P[8] 
RVA/Human-wt/COD/KisB405/2009/G1P[6] 
RVA/Human-wt/COD/KisB501/2009/G1P[6] 
▲ RVA/Human-wt/COD/KisB504/2009/G1 P[6] 
RVA/Human-wt/COD/KisB507/2009/G1P[6] 
RVA/Human-wt/COD/KisB509/2009/G1P[6] 
RVA/Human-wt/COD/KisB512/2009/G1P[6] 
RVA/Human-wt/COD/KisB563/2010/G1P[8] 
ARVA/Human-wt/COD/KisB332/2008/G4P[6] 
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^RVA/Human-wt/COD/KisB603/2008/G4P[6] 
r RVA/Human-wt/COD/KisB609/2008/G12P[6] 
lRVA/Human-wt/COD/KisB610/2008/G12P[6] 
LRVA/Human-wt/COD/KisB323/2008/G12P[6] 
RVA/Human-wt/COD/KisB539/2009/G12P[6] 
RVA/Human-wt/COD/KisB511/2009/G12P[6] 
■RVA/Human-wt/COD/KisB526/2009/G12P[6] 
■RVA/Human-wt/COD/KisB607/2008/G12P[6] 
,RVA/Human-wt/COD/KisB328/2009/G12P[6] 
lRVA/Human-wt/COD/KisB508/2009/G12P[6] 
ARVA/Human-wt/COD/KisB521/2008/G12P[6] 
RVA/Human-wt/COD/KisB534/2009/G12P[6] 
RVA/Human-wt/COD/KisB608/2008/G12P[6] 
RVA/Human-wt/COD/KisB613/2008/G12P[6] 
RVA/Human-wt/COD/KisB331/2009/G9P[8] 
RVA/Human-wt/COD/KisB502/2009/G9P[8] 
RVA/Human-wt/COD/KisB556/2009/G9P[8] 
LRVA/Human-wt/COD/KisB533/2010/G9P[8] 

RVA/Human-wt/GOD/KisB325/2009/G2P[4] 
■RVA/Human-wt/COD/KisB536/2009/G2P[4] 
RVA/Human-wt/COD/KisB527/2009/G2P[4] 
RVA/Human-wt/COD/KisB546/2009/G2P[4] 
RVA/Human-wt/COD/KisB102/2007/G2P[4] 
lRVA/Human-wt/COD/KisB107/2007/G2P[4] 
LRVA/Human-wt/COD/KisB604/2008/G2P[4] 
RVA/Human-wt;COD/KisB523/2009/G2P[4] 
RVA/Human-wt/COD/KisB531/2009/G2P[4] 
RVA/Human-wt/COD/KisB541/2009/G2P[4] 
RVA/Human-wt/COD/KisB542/2009/G2P[4] 
RVA/Human-wt/COD/KisB547/2009/G2P[4] 
RVA/Human-wt/COD/KisB548/2009/G2P[4] 
RVA/Human-wt/COD/KisB558/2009/G2P[4] 
99 r RVA/Human-wt/COD/KisB101/2007/G8P[6] 
LRVA/Human-wt/COD/KisB605/2008/G8P[6] 
RVA/Human-wt/COD/KisB567/2010/G8P[6] 
RVA/Human-wt/COD/KisB561/2010/G8P[6] 
RVA/Human-wt/COD/KisB568/2010/G8P[6] 
RVA/Human-wt/COD/KisB532/2010/G8P[6] 
-RVA/Human-wt/COD/KisB543/2009/G8P[6] 
RVA/Human-wt/COD/KisB545/2009/G8P[6] 
RVA/Human-wt/COD/KisB524/2009/G8P[6] 
RVA/Human-wt/COD/KisB525/2009/G8P[6] 
RVA/Human-wt/COD/KisB529/2010/G8P[6] 
RVA/Human-wt/COD/KisB530/2010/G8P[6] 
RVA/Human-wt/COD/KisB535/2010/G8P[6] 
RVA/Human-wt/COD/KisB538/2009/G8P[6] 
RVA/Human-wt/COD/KisB540/2010/G8P[6] 
RVA/Human-wt/COO/KisB544/2009/G8P[6] 
RVA/Human-wt/COD/KisB549/2010/G8P[6] 
RVA/Human-wt/COD/KisB551/2010/G8P[6] 
RVA/Human-wt/COD/KisB553/2010/G8P[6] 
ARVA/Human-wt/COD/KisB554/2010/G8P[6] 
RVA/Human-wt/COD/KisB555/2010/G8P[6] 
RVA/Human-wt/COD/KisB559/2010/G8P[6] 
RVA/Human-wt/COD/KisB560/2010/G8P[6] 
RVA/Human-wt/COD/KisB562/2010/G8P[6] 
RVA/Human-wt/COD/KisB564/2010/G8P[6] 
ARVA/Human-wt/COD/KisB565/2010/G8P[6] 
RVA/Human-wt/COD/KisB566/2010/G8P[4] 



91 



G1 



G4 



G12 



G9 



G2 
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Figure 5. Maximum likelihood tree based on 363 nucleotides sequence of the VP7 encoding gene of the 99 RVA strains 
characterized in this study. Tree was rooted using RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17] as an outgroup. Bootstrap values (500 replicates) 
above 70 are shown. The strains completely analyzed in this study are marked with a black triangle. 
doi:10.1371/journal.pone.0100953.g005 
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VP6 



97. RVA/Human-wt/ZAF/MRC-DPRU1289/2007/G1P[8] 
» s jRVA/Human-wt/BEL/BE00017/2006/G1P[8] 
jRVA/Human-wt/AUS/CK00083/2008/G1P[8] 
lRVA/Human-wt/BGD/Dhaka12/2003/G12P[6] 
LRVA/Human-wt/DEU/GER172/2008/G12P[6] 
,RVA/Human-wt/KEN/MRC-DPRU2427/2010/G9P[8] 
l-RVA/Human-wt/CMR/MRC-DPRU1417/2009/G9P[8] 
|ARVA/Human-wt/COD/KisB521/2008/G12P[6] 
-' A RVA/Human-wt/COD/KisB504/2009/G1 P[6] 

■ RVA/Human-wt/ZAF/3176WC/2009/G12P[6] 
r RVA/Human-wt/BRA/BA17290/2009/G1 P[6] 
l-RVA/Human-wt/ZWE/MRC-DPRU1723/2009/G9P[8] 
-RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6] 
■ RVA/Human-tc/GBR/ST3/1975/G4P[6] 
LRVA/Human-wt/USA/20 0 8747323/20 08/G1P[8] 

L— RVA/Human-wt/HUN/BP1227/2002/G4P[6] 
RVA/Human-wt/IND/mani-362/2007/G4P[6] 
RVA/Human-wt/ARG/Arg4671/2006/G4P[6] 
RVA/Human-wt/ARG/Arg4605/2006/G4P[6] 
RVA/Human-tc/USA/Wa/1974/G1P[8] 
RVA/Vaccine/USA/Rotarix-RIX4414/1988/G1P1A[8] 
y— RVA/Human-wt/CHN/GX54/2010/G4P[6] 
i — RVA/Human-wt/ECU/EC2184/200x/G11P[6] 
,— RVA/Human-wt/HUN/BP1901/1991/G4P[6] 
l),-RVA/Human-wt/HUN/BP1125/2004/G4P[6] 
M r RVA/Human-wt/HUN/BP1792/2004/G4P[6] 
7 °H r RVA/Human-wt/HUN/BP271/2000/G4P[6] 
, *4-RVA/Human-wt/HUN/BP1490/1994/G4P[6] 
L RVA/Human-wt/HUN/BP1231/2002/G4P[6) 
ARVA/Human-wt/COD/KisB332/2008/G4P[6] 

■ RVA/Pig-tc/USA/Gottfried/1983/G4P[6] 
RVA/Human-wt/NCA/24J/2010/G1P[8] 

RVA/Pig-tc/USA/OSU/1977/G5P9[7] 



-RVA/Human-tc/CHN/R479/2004/G4P[6] 
-RVA/Pig-wt/KOR/K71/2006/G5P[7] 
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100 U 

TT^VA/Human-wt/^ 
^1 — RVA/Human-tc/JPN/AU19/1997/G1P[6] 
-RVA/Human-tc/USA/DS-1/1976/G2P1B[4] 

l RVA/Human-tc/MWI/MW1-860/1999/G8P[6] 

-RVA/Vaccine/USA/RotaTeq-WI79-4/1992/G6P1A[8] 
..rRVA/Human-tc/MWI/QEC29/2005/G8P[6] 
"JlRVA/Human-wt/ZAF/GR10924/1999/G9P[6] 
f-RVA/Human-wt/COD/DRC86/2003/G8P[6] 
\-d.RVA/Human-wt/BEL/B1711/2002/G6P[6] 
4 RVA/Human-wt/BEL/BE1 322/2009/G3P[6] 
••|RVA/Human-wt/USA/06-242/2006/G2P[6] 
lARVA/Human-wt/COD/KisB554/2010/G8P[6] 
„lARVA/Human-wt/COD/KisB565/2010/G8P[6] 
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n A RVA/Human-wt/COD/KisB504/2009/G1 P[6] 
RVA/Human-wt/AUS/CK00083/2008/G1P[8] 
RVA/Human-wt/ZAF/MRC-DPRU1289/2007/G1P[8] 
RVA/Human-wt/BEL/BE00017/2006/G1P[8] 
RVA/Human-wt/ZWE/MRC-DPRU1723/2009/G9P[8] 
RVA/Human-wt/DEU/GER172/2008/G12P[6] 
RVA/Human-wt/CMR/MRC-DPRU1417/2009/G9P[8] 
ARVA/Human-wt/COD/KisB521/2008/G12P[6] 
RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6] 
RVA/Human-wt/KEN/MRC-DPRU2427/2010/G9P[8] 
RVA/Human-wt/ZAF/3176WC/2009/G12P[6l 
RVA/Human-wt/BGD/Dhaka12/2003/G12P[6] 

SRVA/Human-tc/GBR/ST3/1975/G4P[6] 
RVA/Human-wt/USA/20 0 8747323/20 08/G1P[8] 
RVA/Human-tc/USA/Wa/1974/G1P1A[8] 
,„ RVA/Vaccine/USA/Rotarix-RIX4414/1988/G1P1A[8] 
RVA/Human-tc/CHN/R479/2004/G4P[6] 
RVA/Human-wt/CHNK3X54/2010/G4P[6] 
RVA/Human-wt/HUN/BP1901/1991/G4P[6] 
RVA/Human-wt/HUN/BP1490/1994/G4P[6] 
RVA/Human-wt/HUN/BP1792/2004/G4P[6] 
RVA/Human-wt/ARG/Arg4671/2006/G4P[6] 
RVA/Human-wt/ARG/Arg4605/2006/G4P[6] 

ARVA/Human-wt/COD/KisB332/2008/G4P[6] 
RVA/Pig-tc/USA/Gottfried/1983/G4P[6] 
RVA/Pig-tc/USA/OSU/1975/G5P9[7] 
RVA/Pig-wt/KOR/K71/2006/G5P[7] 
RVA/Human-wt/HUN/BP1231/2002/G4P[6] 
RVA/Human-wt/HUN/BP271/2000/G4P[6] 
-RVA/Human-wt/HUN/BP1125/2004/G4P[6] 
RVA/Human-wt/HUN/BP1227/2002/G4P[6] 
RVA/Human-wt/HUN/BP1547/2005/G4P[6] 
RVA/HuSan-tc/JPN/AU19/1997/G1P[6] 
RVA/Vaccine/USA/RotaTeq-WI79-4/1992/G6P1A[8] 
— RVA/Human-tc/USA/DS-1/1976/G2P1B[4] 

ioo,ARVA/Human-wt/COD/KisB554/2010/G8P[6l 
SUn A. RVA/Human-wt/COD/KisB565/2010/G8P 6 
T-RVA/Human-wt/BEL/B1711/2002/G6P[6l 
Kr RVA/Human-wt/COD/DRC86/2003/G8P[6] 
_TlRVA/Human-wt/ZAF/GR10924/1999/G9P[6l 
^1_rRVA/Human-tc/MWI/QEC29/2005/G8P[6] 
, 00 TRVA/Human-tc/MWI/MW1-860/1999/G8P[6] 



R1 



R2 



VP2 




I- A RVA/Human-wt/COD/KisB504/2009/G1 P[6] 
RVA/Human-wt/AUS/CK00083/2008/G1P[8] 
RVA/Human-wt/ZAF/MRC-DPRU1289/2007/G1P[8] 
RVA/Human-wt/BEL/BE00017/2006/G1P[8] 
" RVA/Human-wt/CMR/MRC-DPRU1417/2009/G9P[8] 
RVA/Human-wt/USA/2008747323/2008/G1P[8] 
RVA/Human-wt/BGD/Dhaka12/2003/G12P[6] 
RVA/Human-wt/DEU/GER172/2008/G12P[6] 
• RVA/Human-wt/ZAF/3176WC/2009/G12P[6) 
.-.A RVA/Human-wt/COD/KisB521/2008/G12P[6] 
- B RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6] 
. L RVA/Human-wt/ZWE/MRC-DPRU1723/2009/G9P[8] 
l_ RVA/Human-wt/KEN/MRC-DPRU2427/2010/G9P[8] 
LRVA/Human-tc/GBR/ST3/1975/G4P[6] 
RVA/Human-tc/USA/Wa/1974/G1P1A[8] 
RVA/Vaccine/USA/Rotarix-RIX4414/1988/G1P1A[8] 
RVA/Human-wt/HUN/BP271/2000/G4P[6] 
-RVA/Pig-tc/USA/Gottfried/1983/G4PI6] 
.B,-RVA/Human-wt/HUN/BP1125/2004/G4P[6] 
P— RVA/Human-wt/HUN/BP1231/2002/G4P[6] 

nl RVA/Human-wt/HUN/BP1901/1991/G4P[6] 

| — RVA/Human-wt/HUN/BP1792/2004/G4P[6] 

1 RVA/Human-wt/HUN/BP1547/2005/G4P[6] 

, RVA/Human-wt/ARG/Arg4671/2006/G4P[6] 

|i — RVA/Human-wt/ARG/Arg4605/2006/G4P[6] 

4 RVA/Human-tc/JPN/AU19/1997/G1P[6] 

, ARVA/Human-wt/COD/KisB332/2008/G4P[6] 

L.RV - 
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RVA/Pig-tc/USA/OSU/1977/G5P9[7] 
s ,-RVA/Pig-wt/KOR/K71/2006/G5P[7] 
RVA/Human-wt/HUN/BP1227/2002/G4P[6] 
RVA/Human-tc/CHN/R479/2004/G4P[6] 
-RVA/Human-wt/CHN/GX54/2010/G4P[6] 

RVA/Vaccine/USA/RotaTeq-WI79-4/1992/G6P1A[8] 

RVA/Human-tc/USA/DS-1/1976/G2P1B[4] 

.0j-RVA/Human-tc/MWI/QEC29/2005/G8P[6] 
jtRVA/Human-wt/ZAF/GR10924/1999/G9P[6] 
L.RVA/Human-wt/COD/DRC86/2003/G8P[6] 
RVA/Human-tc/MWI/MW1-850/1999/G8P[6] C2 
RVA/Human-wt/BEL/B1711/2002/G6P[6] 
RVA/Human-wt/BEL/BE1322/2009/G3P[6] 
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C1 



Figure 6. Maximum likelihood trees inferred from the nucleotide sequences of RVA VP6, VP1 and VP2 genes. The gene sequence 
lengths compared were: VP6 (1064 bp), VP1 (3168 bp), VP2 (2534 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values < 
70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17]) for each 
gene was included. Completely sequenced RVA strains from Kisangani are indicated with black triangles. 
doi:1 0.1 371 /journal.pone.01 00953.g006 
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5 RVA/Vaccine/USA/Rotarix-RIX4414/1988/G1P1A[8] 
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-RVA/Human-tc/JPN/AU19/1997/G1P[6] 
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RVA/Human-wt/ZMB/MRC-DPRU3491/2009/G12P[6] 
ARVA/Human-wt/COD/KisB521/2008/G12P[6] 
jARVA/Human-wt/COD/KisB504/2009/G1P[6l 
RVA/Human-wt/DEU/GER172/2008/G12P[6] 
RVA/Human-wt/ZAF/3176WC/2009/G12P[6] 
RVA/Human-wt/BRA/BA17290/2009/G1P[6] 
RVA/Human-wt/USA/20 08747323/20 08/G1P[8] 
RVA/Human-wt/KEN/MRC-DPRU2427/2010/G9P[8] 
RVA/Human-tc/GBR/ST3/1975/G4P[6] 

„, — RVA/Human-wt/BRA/NB150/1997/G1G4P[6l 

J lRVA/Vaccine/USA/Rotarix-RIX4414/1988/G1P1A[8] 

I RVA/Human-tc/USA/Wa/1 974/G1 P1A[8] 

ARVA/Human-wt/COD/KisB332/2008/G4P[6] 
RVA/Human-tc/USA/WI61/1983/G9P1A[8] 
■ RVA/Pig-tc/USA/OSU/1977/G5P9[7] 
,„ RVA/Pig-wt/KOR/K71/2006/G5P[7] 
RVA/Human-wt/HUN/BP1231/2002/G4P[6] 
RVA/Human-tc/CHN/R479/2004/G4P[6] 
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RVA/Human-wt/HUN/BP1227/2002/G4P[6] 
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Figure 7. Maximum likelihood trees inferred from the nucleotide sequences of RVA VP3, NSP1 and NSP2 genes. The gene sequence 
lengths compared were: VP3 (2442 bp), NSP1 (1 1 39 bp), NSP2 (938 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values 
<70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17] for VP3 
and NSP2 and RVA/Horse-wt/ARG/E403/2006/G14P[12] for NSP1) was included. Completely sequenced RVA strains from Kisangani are indicated with 
black triangles. 

doi:10.1371/journal.pone.0100953.g007 
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Figure 8. Maximum likelihood trees inferred from the nucleotide sequences of RVA NSP3, NSP4 and NSP5 genes. The gene sequence 
lengths compared were: NSP3 (841 bp), NSP4 (501 bp), NSP5 (542 bp). The bootstrap values (500 replicates) are shown at the branch nodes (values < 
70% not shown). The scale bar is proportional to the genetic distance. An outgroup sequence (RVA/Pigeon-tc/JPN/PO-13/1983/G18P[17]) for each 
gene was included. Completely sequenced RVA strains from Kisangani are indicated with black triangles. 
doi:10.1371/journal.pone.0100953.g008 
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investigating 6 Brazilian G1P[6] RVA strains, isolated from 
vaccinated children in 2009 and 2010 [42]. The G12-P[6]-I1-R1- 
C1-M1-A1-N1-T1-E1-H1 configuration has been previously 
reported in two African countries: South-Africa (RVA/Human- 
wt/ZAF/3176WC/2009/G12P[6]) and Zambia (RVA/Human- 
wt/ZMB/MRC-DPRU3491/2009/G12P[6]) [35]. Both G8P[6] 
strains possessed a conserved genomic configuration: G8-P[6]T2- 
R2-C2-M2-A2-N2-T2-E2-H2. This constellation was exactly the 
same as that of RVA/Human-wt/COD/DRC86/2003/G8P[6], 
also isolated in the COD, and RVA/Human-tc/MWI/QEC29/ 
2005/G8P[6], isolated in Malawi [11,13]. 

The G4P[6] strain was found to possess the, for humans 
unusual, G4-P[6]-I1-R1-C 1-M1-A1-N1-T7-E1-H1 genotype con- 
stellation. Although unusual in humans, a recent study describing 
G4P[6] strains in Hungarian children, reported three strains with 
the exact same genotype constellation as KisB332 [43]. The 
genotype constellation I1/I5-R1-C1-M1-A1/A8-N1-T1/T7-E1- 
H 1 is typically found in pigs with a variety of different G- and P- 
genotypes [9,44-46] . A number of recent studies have identified 
human G4P[6] strains detected in different parts of the world, 
which are believed to be examples of direct interspecies 
transmissions, or may have undergone reassortments with human 
RVA strains [43,47-51]. 

Phylogenetic analyses. Despite the different G-genotype, a 
very close genetic relationship exists between the G1P[6] 
(KisB504) and G12P[6] (KisB521) strain within their segments 
VP4 (99.6%), VP6 (99.9%), NSP1 (99.6%), NSP4 (99.8%) and 
NSP5 (99.8%). For the remaining genome segments encoding 
VP1, VP2, VP3, NSP2 and NSP3 a slighdy lower nucleotide 
sequence similarity was detected, ranging from 94.4% (for NSP3) 
to 98.2% (for NSP2). The genetic similarity between the VP 1-4, 
VP6, VP7, NSP1-3 and NSP5 of RVA strains KisB521, KisB504 
and other recently isolated human Wa-like RVA strains from all 
over the world was very high (figure 4-8). The majority of gene 
segments of strain KisB521 (VP7, VP4, VP1-3, NSP1, NSP4) and 
five gene segments of strain KisB504 (VP4 and NSP1-4) are closely 
related to a recently isolated strain from Zambia (RVA/Human- 
wt/ZMB/MRC-DPRU3491/2009/G12P[6]). 

Both G8P[6] strains clustered very closely together in the 
phylogenetic trees of all gene segments, with a nucleotide sequence 
similarity ranging between 99.5% and 100% (figure 4-8). These 
results indicate a clonal origin of both strains, which were isolated 
within a 2-month interval, from patients admitted to the same 
hospital. For VP7, VP3, NSP3 and NSP5 these strains shared the 
highest nucleotide similarity with RVA/Human-wt/COD/ 
DRC86/2003/G8P[6]. Human RVA strain DRC86 is the only 
completely characterized Congolese G8P[6] strain available in 
Genbank until now [13]. Maximal nucleotide sequence identity of 
KisB554 and KisB565 was found with RVA strain RVA/Human- 
wt/USA/06-242/2006/G2P[6] (VP4, VP6, VP2, NSP1 and 
NSP2) and RVA/Human-wt/BEL/B 1711/ 2002/ G6P [6] (VP1). 
Strain 06-242 was isolated in an American child in 2006 and was 
the first G2P[6] strain characterized with a DS-l-like genotype 
constellation [52]. Strain B1711 also possessed a complete DS-l- 
like genotype constellation and was isolated from a 13-month-old 
child admitted to the university hospital in Leuven, Belgium, with 
severe gastroenteritis after returning from a trip to Mali [53,54]. 

Overall, most gene segments of the G8P[6], G1P[6] and 
G12P[6] strains appeared to cluster closely to gene segments of 
recently isolated human RVA strains from all around the world, 
with the exception of NSP4. This gene segment showed more 
similarity with porcine/human porcine-like (for the G1P[6], 
G12P[6] and G4P[6] strains) or bovine/human bovine-like (for 
the G8P[6] strains) RVAs (figure 8). More specific, the three 



Congolese strains (and the Zimbabwean strain RVA/Human-wt/ 
ZMB/MRC-DPRU3491/2009/G12P[6]) belonging to genotype 
El clustered together with RVA strains RVA/Human-wt/BRA/ 
NB150/1997/G1G4P[6] and RVA/Human-wt/ECU/EC2184/ 
200x/Gl 1P[6], both previously classified as porcine-human 
reassortants [51,55]. Furthermore, these six strains also felt into 
a larger cluster containing porcine RVA strains, or strains believed 
to be of porcine origin. On the other hand, the NSP4 encoding 
gene segments of the G8P[6] strains belonged to an E2 genotype 
cluster containing both animal and human P[14] RVA strains, 
which are believed to be the result of direct interspecies 
transmissions from sheep or other members of the mammalian 
order Artiodactyla [6] . 

For the G4P[6] strain (KisB332) it appeared that all gene 
segments were more closely related to porcine RVA strains or 
human porcine-like RVA strains than to typical human RVA 
strains. A possible exception is VP7, because it is difficult to 
determine the exact origin of the VP7 gene, since there are no 
porcine G4 RVA sequences from Africa known today. In general, 
the gene segments of KisB332 showed lower nucleotide sequence 
identities to known contemporary and historic human or animal 
RVAs, especially in the VP6, VP1, VP2 and NSP1 genes (figures 6 
and 7). Although these gene segments of KisB332 were only 
relatively distantly related to known RVA strains, they showed the 
highest relatedness to RVA/Human-wt/NCA/24J/2010/GlP[8] 
(for VP6), RVA/Pig-tc/USA/OSU/1977/G5P9[7] (VP1), RVA/ 
Pig-wt/KOR/K71/2006/G5P7 (VP2) and human reference 
strain RVA/Human-tc/USA/WI61/1983/G9PlA[8] (NSP1) 
[56,57]. The VP3 of KisB332 was most closely related to porcine 
RVA strain OSU (figure 7). The VP4, NSP2, NSP3, NSP4 and 
NSP5 genes clustered together with putative porcine interspecies 
transmitted RVA strains or human-porcine reassortant strains 
(BE2001, Arg4671, NB150, EC2184) [47,51,55,58]. Interestingly, 
the NSP3 gene segment of KisB332 belongs to the rarely described 
T7 genotype (figure 8). The T7 genotype was first detected in a 
bovine RVA strain (RVA/Cow-tc/GBR/UK/1973/G6P7[5]) 
and was more recendy described in a bovine-like human strain 
(RVA/Human-xx/IND/Mani-265/2007/G10P[6]) [48,59]. 
However, an increasing number of studies have reported the T7 
genotype in porcine-like human strains (R479, BE2001, BP 1901, 
BP1490, BP271, BP1792, BP1125, BP1547 and Arg4605) and 
porcine RVA strains (CMP45, AB82, F6-4, F8-4-A and TM-a) 
derived from different geographical locations [43,46,47,50,58]. 
The detection of the T7 genotype in strain RVA/Human-wt/ 
COD/KisB332/2008/G4P[6], a strain with a human porcine-like 
gene constellation, makes it more likely that T7 genotypes have 
their origins in pigs as was previously suggested in other studies 
[46,58]. 

Discussion 

Data shown in figure 2, suggest three periods of increased RVA 
detection rates over the course of this study. However, these 
fluctuations can be partially explained by sampling bias inherent to 
the circumstances in which these samples have been collected. 
Kisangani is known for its complex political and economical 
climate, making it difficult to collect sufficient samples over a large 
period of time without any period of lower or no sampling. 
Although the data available in this study do have their limitations, 
the available data confirm previous studies that reported year- 
round RVA infections in tropical settings, with ups and downs but 
no clear period with minimal or no RVA hospitalizations [60,61]. 

The six most common G/P-genotype combinations in countries 
with a more moderate climate and by extension worldwide are 
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G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8]. In this 
study a high proportion of P[6] RVA strains was observed together 
with a large number of unusual strains (G8P[6], G12P[6], G1P[6], 
G4P[6] and G8P[4]). These unusual genotype combinations may 
represent potential reassortants between human and animal RVA 
strains as a result of zoonotic infection. The only way to investigate 
this is to perform phylogenetic analysis of complete genomes. 
Because of the high number of P[6] RVA strains detected in this 
study and the fact that the P[6] genotype is much less prevalent in 
high-income countries this study focused on the four different P[6] 
genotype combinations detected in Kisangani between May 2007 
and December 2010. The majority of the gene segments from the 
five P[6] strains characterized in this study appeared to have a 
genetic background most closely related to human RVA strains 
circulating all over the world. However, all five strains possessed at 
least one gene segment (NSP4) that was most closely related to 
animal RVAs. One strain, RVA/Human-wt/COD/KisB332/ 
2008/G4P[6], possessed a majority of gene segments (if not all) 
most closely related to porcine, or human porcine-like RVA 
strains. For most genome segments KisB332 had still a consider- 
able genetic distance to known porcine RVA strains. This can 
partly be explained by sampling bias since surveillance studies for 
animal rotaviruses are much less common than human rotavirus 
surveillance studies. In addition, porcine RVA sequences obtained 
from African pigs are not available in Genbank thus far. The lack 
of availability of complete genome sequences from porcine RVA 
strains in general and from African pig RVAs in particular makes 
it difficult to differentiate between the porcine or human origin of 
certain gene segments. 

The child from which strain KisB332 was isolated lived in close 
proximity of a pig farm, providing a plausible route of 
transmission. Unfortunately it was not possible to take samples 
from the pigs living in close contact with the child. These findings 
suggest that KisB332 is of porcine origin and that this strain, after 
an interspecies transmission and potential reassortment with 
human RVAs, was able to cause gastro-enteritis in a human 
child. Surprisingly, a second G4P[6] strain (KisB603) was found in 
this study with very similar VP4 and VP7 gene sequences as those 
of KisB332. In theory, both strains can be the result of 
independent interspecies transmission events but a more plausible 
theory would be that the G4P[6] strain was able to spread from 
one human to another, although this hypothesis needs further 
investigation. 

Although interspecies transmission events mostly result in so- 
called death-end infections, interspecies transmissions from pigs to 
humans seem to occur rather frequently, probably due to the fact 
that porcine RVA strains share a similar genetic background with 
human Wa-like RVAs [62]. Overall, these results underline the 
importance of complete genetic characterization of RVA strains 

References 

1. Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Duquc J, et al. (2012) the 
WHO-eoordinated Global Rotavirus Surveillanee Network. 2008 estimate of 
worldwide rotavirus-assoeiated mortality in children younger than 5 years before 
the introduction of universal rotavirus vaccination programmes: a systematic 
review and meta-analysis. Lancet Infect Dis 12: 136-141. 

2. Schcri M, Nemarude L, Peenzc I, Nctshiihcihc L, Nyaga MM, ct al. (2014) 
Update of rotavirus strains circulating in Africa from 2007 through 20 1 1 . Pediatr 
Infect Dis J 33 Suppl 1: S76-84. 

3. Jerc KC, Mlera L, O'Neill HG, Peenzc I, van Dijk AA (2012) Whole genome 
sequence analyses of three African bovine rotaviruses reveal that they emerged 
through multiple reassortment events between rotaviruses from different 
mammalian species. Vet Microbiol 159: 245-250. 

4. Alkan F, Ozkul A, Oguzoglu TC, Timurkan MO, Caliskan E, et al. (2010) 
Distribution of G (VP7) and P (VP4) genotypes of group A bovine rotaviruses 
from Turkish calves with diarrhea, 1997-2008. Vet Microbiol 141: 231-237. 



and indicate that reassortments and interspecies transmission 
between human RVAs and porcine RVA strains occur frequentiy 
in the COD. Interspecies transmission events can have a potential 
negative influence on the efficacy of RVA vaccines, highlighting 
the need for continued surveillance of RVA diversity in the COD. 

Based on the data from this study, RVA vaccines will be 
challenged with a wider variety of different G- and P- genotypes in 
the COD compared to those found in high-income countries. Of 
particular importance for the potential implementation of the 
currently available RVA vaccines in the COD is the predomi- 
nance of P[6] detected in this study (52%) and the high proportion 
of strains with the G8 VP7 specificity (27%). Since these genotypes 
are not included in the currendy available RVA vaccines and the 
genotype specific vaccine efficacy against these genotypes is not yet 
known, the success of vaccination programs could be lower. 
Nevertheless, implementation of rotavirus vaccines in the COD 
would be highly recommended to prevent rotavirus related 
morbidity and mortality. 

Supporting Information 

Figure SI. Phylogenetic dendrogram based on 595 
nucleotides of the VP4 encoding gene of the 12 repre- 
sentative RVA strains characterized in the study (black 
triangle) and RVA strains isolated all over the world. 

Tree was rooted using RVA/Human-tc/JPN/AU-l/1982/G3P9] 
as an outgroup. Bootstrap values (500 replicates) above 70 are 
shown. 
(TIF) 

Figure S2 Phylogenetic dendrogram based on 503 
nucleotides of the VP7 encoding gene of the 12 repre- 
sentative RVA strains characterized in this study (black 
triangle) and RVA strains isolated all over the world. 

Tree was rooted using RVA/Pigeon-tc/JPN/PO-13/1983/ 
G18P[17] as an outgroup. Bootstrap values (500 replicates) above 
70 are shown. 
(TIF) 

Table SI Primers used to amplify the VP1, VP2, VP3, 
VP4, VP6, VP7, NSP1, NSP2, NSP3, NSP4 and NSP5 gene 
segments described in this study. 

(DOC) 

Author Contributions 

Conceived and designed the experiments: EH BBL CVG JJ DN JM. 
Performed the experiments: EH BBL MZ SS SDC. Analyzed the data: EH 
MZ NGNJM. Contributed reagents/materials/analysis tools: MVR CVG 
DN. Wrote the paper: EH JM. 



5. Midgley SE, Banyai K, Buesa J, Halaihcl N, Hjulsager CK, et al. (2012) 
Diversity and zoonotic potential of rotaviruses in swine and cattle across Europe. 
Vet Microbiol 156: 238-245. 

6. Matthijnssens J, Potgietcr CA, Ciarlct M, Parreno V, Martclla V, ct al. (2009) 
Are human P[14J rotavirus strains the result of interspecies transmissions from 
sheep or other ungulates that belong to the mammalian order Artiodactyla? 
J Virol 83: 2917-2929. 

7. Monini M, Cappuccini F, Battista P, Falcone E, Lavazza A, ct al. (2008) 
Molecular characterization of bovine rotavirus strains circulating in northern 
Italy, 2003-2005. Vet Microbiol 129: 384-389. 

8. Cunliffe NA, Gentsch JR, Kirkwood CD, Gondwe JS, Dove W, ct al. (2000) 
Molecular and serologic characterization of novel serotype G8 human rotavirus 
strains detected in Blantyre, Malawi. Virology 274: 309—320. 

9. Papp H, Laszlo B, Jakab F, Gancsh B, Dc Grazia S, ct al. (2013) Review of 
group A rotavirus strains reported in swine and cattle. Vet Microbiol 165: 190- 
199. 



PLOS ONE | www.plosone.org 



12 



June 2014 | Volume 9 | Issue 6 | e100953 



Rotavirus Surveillance in Kisangani 



10. Martclla V, Banyai K, Ciarlct M, Iturriza-Gomara M, Lorusso E, et al. (2006) 
Relationships among porcine and human P[6J rotaviruses: evidence that the 
different human P[6] lineages have originated from multiple interspecies 
transmission events. Virology 344: 509-519. 

11. Nakagomi T, Doan YH, Dove W, Ngwira B, Iturriza-Gomara M, et al. (2013) 
G8 rotaviruses with conserved genotype constellations detected in Malawi over 
10 years (1997-2007) display frequent gene reassortmcnt among strains co- 
circulating in humans. J Gen Virol 94: 1273—1295. 

12. Esona MD, Geyer A, Page N, Trabelsi A, Fodha I, et al. (2009) Genomic 
characterization of human rotavirus G8 strains from the African rotavirus 
network: relationship to animal rotaviruses. J Med Virol 81: 937—951. 

13. Matthijnsscns J, Rahman M, Yang X, Delbeke T, Arijs I, et al. (2006) G8 
rotavirus strains isolated in the Democratic Republic of Congo belong to the DS- 
1-like genogroup. J Clin Microbiol 44: 1801-1809. 

14. Kabuc JP, Pecnze I, de Beer M, Esona MD, Lunfungula C, et al. (2010) 
Characterization of human rotavirus recovered from children with acute 
diarrhea in Kinshasa, Democratic Republic Of Congo. J Infect Dis 202 Suppl: 
S193-197. 

15. Nielsen NM, Eugcn-Olsen J, Aaby P, Molbak K, Rodrigues A, ct al. (2005) 
Characterisation of rotavirus strains among hospitalised and non-hospitalised 
children in Guinea-Bissau, 2002 A high frequency of mixed infections with 
serotype G8. J Clin Virol 34: 13-21. 

16. Nokcs DJ, Peenzc I, Netshifhefhe L, Abwao J, De Beer MC, et al. (2010) 
Rotavirus genetic diversity, disease association, and temporal change in 
hospitalized rural Kenyan children. J Infect Dis 202 Suppl: S180— 186. 

17. Cunliffe NA, Ngwira BM, Dove W, Thindwa BD, Turner AM, et al. (2010) 
Epidemiology of rotavirus infection in children in Blantyre, Malawi, 1997—2007. 
J Infect Dis 202 Suppl: SI 68- 174. 

18. Aminu M, Page NA, Ahmad AA, UmohJU, Dewar J, et al. (2010) Diversity of 
rotavirus VP7 and VP4 genotypes in Northwestern Nigeria. J Infect Dis 202 
Suppl: S198-204. 

19. Jerc KC, Sawyerr T, Schcri LM, Pecnze I, Page NA, et al. (201 1) A first report 
on the characterization of rotavirus strains in Sierra Leone. J Med Virol 83: 540- 
550. 

20. Binka EN, Anto EK, Oduro AR, Awini EA, Nazzar AK, et al. (2003) Incidence 
and risk factors of paediatric rotavirus diarrhoea in northern Ghana. Trop Med 
Int Health 8: 840-846. 

21. Yassin MA, Kirby A, Mengistu AA, Arbide I, Dove W, et al. (2012) Unusual 
norovirus and rotavirus genotypes in Ethiopia. Paediatr Int Child Health 32: 51- 
55. 

22. Ndze VN, Papp H, Achidi EA, Gonsu KH, Laszlo B, et al. (2013) One year 
survey of human rotavirus strains suggests the emergence of genotype G12 in 
Cameroon. J Med Virol 85: 1485-1490. 

23. Nordgren J, Bonkoungou IJ, Nitiema LW, Sharma S, Ouermi D, et al. (2012) 
Rotavirus in diarrheal children in rural Burkina Easo: high prevalence of 
genotype G6P[6]. Infect Genet Evol 12: 1892-1898. 

24. Ruiz-Palacios GM, Perez-Schael I, Velazquez ER, Abate H, Breuer T, et al. 
(2006) Safety and efficacy of an attenuated vaccine against severe rotavirus 
gastroenteritis. N Engl J Med 354: 1 1-22. 

25. Vesikari T, Matson DO, Dennchy P, Van Damme P, Santosham M, ct al. 
(2006) Safety and efficacy of a pcntavalcnt human -bovine (WC3) reassortant 
rotavirus vaccine. N Engl J Med 354: 23—33. 

26. Armah GE, Sow SO, Breiman RF, Dallas MJ, Tapia MD, et al. (2010) Efficacy 
of pentavalent rotavirus vaccine against severe rotavirus gastroenteritis in infants 
in developing countries in sub-Saharan Africa: a randomised, double-blind, 
placebo-controlled trial. Lancet 376: 606—614. 

27. Madhi SA, Cunliffe NA, Steele D, Witte D, Kirsten M, et al. (2010) Effect of 
human rotavirus vaccine on severe diarrhea in African infants. N Engl J Med 
362: 289-298. 

28. Sow SO, Tapia M, Haidara FC, Ciarlet M, Diallo F, et al. (2012) Efficacy of the 
oral pcntavalcnt rotavirus vaccine in Mali. Vaccine 30 Suppl 1: A71-78. 

29. Johns Hopkins Bloomberg School of Public Health International Vaccine Access 
Center. (2014) VIMS Report: Global Vaccine Introduction. Database: VIMS 
http:/ / www.jhsph.edu/ research/ centers-and-institutcs/ivac/vims/IVAC- 
VIMS-Report-2014-Mar.pdf. Accessed 4June 2014. 

30. PATH (2014) Rotavirus vaccines access and delivery. Country introductions of 
rotavirus vaccines. Available: http://sites.path.org/rotavirusvaccine/rotavirus- 
advocacy-and-communieations-toolkit/ country-introduction-maps-and-list/. 
Accessed 4 June 2014. 

31. Moon SS, Tate JE, Ray P, Dennehy PH, Archary D, et al. (2013) Differential 
Profiles and Inhibitory Effect on Rotavirus Vaccines of Non-Antibody 
Components in Breast Milk from Mothers in Developing and Developed 
Countries. Pediatr Infect Dis J. 

32. Serazin AC, Shackclton LA, Wilson C, Bhan MK (2010) Improving the 
performance of enteric vaccines in the developing world. Nat Immunol 1 1: 769— 
773. 

33. Hu L, Crawford SE, Czako R, Cortes-Penfield NW, Smith DE, et al. (2012) Cell 
attachment protein VP8* of a human rotavirus specifically interacts with A-type 
histo-blood group antigen. Nature 485: 256-259. 

34. Nyaga MM, Jcre KC, Peenze I, Mlcra L, van Dijk AA, ct al. (2013) Sequence 
analysis of the whole genomes of live African human G9 rotavirus strains. Infect 
Genet Evol 16: 62-77. 

35. Jerc KC, Mlcra L, O'Neill HG, Potgieter AC, Page NA, et al. (2011) Whole 
genome analyses of African G2, G8, G9, and G12 rotavirus strains using 



sequence-independent amplification and 454(R) pyroscquencing. J Med Virol 
83: 2018-2042. 

36. Esona MD, Mijatovic-Rustempasic S, Foytich K, Roy S, Banyai K, et al. (2013) 
Human G9P[8] rotavirus strains circulating in Cameroon, 1999-2000: Genetic 
relationships with other G9 strains and detection of a new G9 subtype. Infect 
Genet Evol 18C: 315-324. 

37. Ghosh S, Gatheru Z, NyangaoJ, Adachi N, Urushibara N, et al. (2011) Full 
genomic analysis of a G8P[1J rotavirus strain isolated from an asymptomatic 
infant in Kenya provides evidence for an artiodactyl-to-human interspecies 
transmission event. J Med Virol 83: 367—376. 

38. Jerc KC, Mlera L, Page NA, van Dijk AA, O'Neill HG (201 1) Whole genome 
analysis of multiple rotavirus strains from a single stool specimen using sequence- 
independent amplification and 454(R) pyroscquencing reveals evidence of 
intergenotype genome segment recombination. Infect Genet Evol 11: 2072- 
2082. 

39. Rahman M, Goegebuer T, De Leener K, Maes P, Matthijnsscns J, ct al. (2004) 
Chromatography paper strip method for collection, transportation, and storage 
of rotavirus RNA in stool samples. J Clin Microbiol 42: 1605—1608. 

40. Matthijnssens J, Rahman M, Martclla V, Xuelei Y, De Vos S, ct al. (2006) Full 
genomic analysis of human rotavirus strain B4106 and lapinc rotavirus strain 
30/96 provides evidence for interspecies transmission. J Virol 80: 3801-3810. 

41. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, evolution- 
ary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731-2739. 

42. Gomez MM, da Silva MF, Zeller M, Heylcn E, Matthijnssens J, ct al. (2013) 
Phylogenetie analysis of GlP[6j group A rotavirus strains detected in Northeast 
Brazilian children fully vaccinated with Rotarix. Infect Genet Evol 19: 395-402. 

43. Papp H, Borzak R, Farkas S, Kisfali P, Lengyel G, ct al. (2013) Zoonotic 
transmission of reassortant porcine G4P[6j rotaviruses in Hungarian pediatric 
patients identified sporadically over a 15 year period. Infect Genet Evol 19: 71- 
80. 

44. Kim HH, Matthijnssens J, Kim HJ, Kwon HJ, ParkJG, et al. (2012) Full-length 
genomic analysis of porcine G9P[23] and G9P[7J rotavirus strains isolated from 
pigs with diarrhea in South Korea. Infect Genet Evol 12: 1427-1435. 

45. Okitsu S, Khamrin P, Thongprachum A, Kongkacw A, Maneekarn N, et al. 
(2013) Whole-genomic analysis of G3P[23], G9P[23J and G3P[13] rotavirus 
strains isolated from piglets with diarrhea in Thailand, 2006-2008. Infect Genet 
Evol 18: 74-86. 

46. Martcl-Paradis O, Laurin MA, Martella V, SohalJS, L'Homme Y (2013) Full- 
length genome analysis of G2, G9 and Gil porcine group A rotaviruses. Vet 
Microbiol 162: 94-102. 

47. DegiuscppcJI, Beltramino JC, Millan A, StupkaJA, Parra GI (2013) Complete 
genome analyses of G4P[6J rotavirus detected in Argentinean children with 
diarrhoea provides evidence of interspecies transmission from swine. Clin 
Microbiol Infect 19: E367-371. 

48. Mukherjee A, Ghosh S, Bagchi P, Dutta D, Chattopadhyay S, ct al. (201 1) Full 
genomic analyses of human rotavirus G4P[4], G4P[6], G9P[19] and G10P[6J 
strains from North-eastern India: evidence for interspecies transmission and 
complex reassortmcnt events. Clin Microbiol Infect 17: 1343-1346. 

49. Dong HJ, Qjan Y, Huang T, Zhu RN, Zhao LQ, et al. (2013) Identification of 
circulating porcine-human reassortant G4P[6] rotavirus from children with 
acute diarrhea in China by whole genome analyses. Infect Genet Evol 20: 155— 
162. 

50. Wang YH, Kobayashi N, Nagashima S, Zhou X, Ghosh S, et al. (2010) Full 
genomic analysis of a porcine-bovine reassortant G4P[6] rotavirus strain R479 
isolated from an infant in China. J Med Virol 82: 1094-1102. 

51. Maestri RP, Kaiano JH, Neri DL, Soares Lda S, Guerra Sde F, et al. (2012) 
Phylogenetie analysis of probable non-human genes of group A rotaviruses 
isolated from children with acute gastroenteritis in Belcm, Brazil. J Med Virol 
84: 1993-2002. 

52. Hcylen E, Zeller M, Ciarlet M, De Coster S, Van Ranst M, et al. (2013) 
Complete genetic characterization of human G2P[6] and G3P[6J rotavirus 
strains. Infect Genet Evol 13: 27—35. 

53. Rahman M, De Lcencr K, Goegebuer T, Wollants E, Van dcr Donck I, et al. 
(2003) Genetic characterization of a novel, naturally occurring recombinant 
human G6P[6J rotavirus. J Clin Microbiol 41: 2088-2095. 

54. Matthijnsscns J, Rahman M, Van Ranst M (2008) Two out of the 1 1 genes of an 
unusual human G6P[6J rotavirus isolate arc of bovine origin. J Gen Virol 89: 
2630-2635. 

55. Banyai K, Esona MD, Kcrin TK, Hull JJ, Mijatovic S, et al. (2009) Molecular 
characterization of a rare, human-porcine reassortant rotavirus strain, Gl 1P[6J, 
from Ecuador. Arch Virol 154: 1823-1829. 

56. Bucardo F, Rippinger CM, Svcnsson L, PattonJT (2012) Vaccine-derived NSP2 
segment in rotaviruses from vaccinated children with gastroenteritis in 
Nicaragua. Infect Genet Evol 12: 1282-1294. 

57. Clark HF, Hoshino Y, Bell LM, GroffJ, Hess G, et al. (1987) Rotavirus isolate 
WI61 representing a presumptive new human serotype. J Clin Microbiol 25: 
1757-1762. 

58. Zeller M, Hcylen E, Dc Coster S, Van Ranst M, Matthijnssens J (2012) Full 
genome characterization of a porcine-like human G9P[6] rotavirus strain 
isolated from an infant in Belgium. Infect Genet Evol 12: 1492-1500. 

59. Bridger JC, Woode GN (1975) Neonatal calf diarrhoea: identification of a 
rcovirus-likc (rotavirus) agent in faeces by immunofluorescence and immune 
electron microscopy. Br Vet J 131: 528-535. 



PLOS ONE | www.plosone.org 



13 



June 2014 | Volume 9 | Issue 6 | e100953 



Rotavirus Surveillance in Kisangani 



60. Cook SM, Glass RI, LcBaron GW, Ho MS (1990) Global seasonality of rotavirus 
infections. Bull World Health Organ 68: 171-177. 

6 1 . Levy K, Hubbard AE, Eisenberg JN (2009) Seasonality of rotavirus disease in 
the tropics: a systematic review and meta-analysis. Int J Epidemiol 38: 1487- 
1496. 



62. Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Dclbckc T, ct al. (2008) Full 
genome-based classification of rotaviruses reveals a common origin between 
human Wa-Like and porcine rotavirus strains and human DS-l-like and bovine 
rotavirus strains. J Virol 82: 3204-3219. 



PLOS ONE | www.plosone.org 



14 



June 2014 | Volume 9 | Issue 6 | e100953 



